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Abstract — Electric  potentials  across  the  auro-al  zones  and  polar  cap  have  been  determined  by  an  integration 
of  electric  fields  measured  on  the  Dynamics  explorer  2  satellite.  Potentials  from  527  cases,  divided  between 
winter  and  summer  seasons  in  the  northern  hemisphere,  are  compared  with  simultaneous  aurora]  electrojet 
indices  which  are  measured  with  ground-based  magnetic  observatories.  The  dusk  and  dawn  auroral  potentials 
and  the  total  polar  cap  potential  are  compared  with  the  AU,  AL,  and  AE  indices  respectively.  The  AE 
■ndices  are  found  to  be  higher  when  it  is  summer  in  the  northern  hemisphere,  due  to  the  higher  ionospheric 
conductivity.  The  electrojet  currents  are  roughly  proportional  to  the  potentials,  but  the  AE  indices  have  a 
range  of  possible  values  for  a  given  potential  due  to  variations  in  the  width  and  distribution  of  the  electrojets, 
the  ionospheric  Hall  conductivity,  and  the  neutral  wind  velocity.  The  data  have  been  fit  to  linear  and  non¬ 
linear  equations,  which  may  be  used  to  estimate  the  potentials  from  the  AE  indices.  Although  the  non-linear 
equations  give  slightly  better  fits  than  a  straight  line,  the  differences  in  the  errors  are  not  significant.  Where 
possible  the  data  have  also  been  grouped  according  to  the  phase  of  magnetospheric  substorms.  These  results 
indicate  that  substorms  occur  only  when  the  polar  cap  potential  exceeds  a  threshold  of  approximately  60  kV. 


1.  INTRODUCTION 

In  this  paper  we  will  show  the  results  of  a  compari¬ 
son  of  auroral  and  polar  cap  electric  potentials  and 
the  Auroral  Electrojet  (AE)  indices.  The  electric 
potential  measurements  are  from  the  Vector  Electric 
Field  Instrument  (VEFI)  on  the  Dynamics  Explorer- 
2  (DE-2)  satellite.  Potentials  measured  across  the 
dawn  and  dusk  auroral  zones  and  the  polar  cap  are 
compared  with  the  AL,  AU,  and  AE  indices  at  the 
time  of  the  satellite  measurements. 

There  are  several  reasons  why  this  comparison  be¬ 
tween  potentials  and  AE  indices  can  be  useful.  First 
of  all,  knowing  how  the  potentials  and  electrojet  cur¬ 
rents  relate  to  each  other  will  promote  an  under¬ 
standing  of  the  electrodynamics  of  the  aurora  un¬ 
der  both  steady  state  and  and  substorm  conditions. 
Secondly,  although  there  have  been  studies  of  the 
relations  between  the  interplanetary  magnetic  field 
(IMF)  and  the  AE  index  and  between  the  IMF  and 
polar  cap  potential  (in  separate  studies),  there  has 
never  been  an  in-depth  comparison  between  AE  and 
potentials.  With  a  functional  relationship  which  is 


derived  from  simultaneous  measurements  of  the  po¬ 
tential  by  a  satellite  and  the  magnetic  perturbations 
on  the  ground,  it  will  be  possible  to  estimate  the 
potentials  from  the  magnetic  activity  indices.  This 
could  be  valuable  to  investigators  who  make  use  of 
the  various  electric  field  models;  it  would  be  useful 
to  know  if  the  total  electric  potential  in  a  model  is  an 
accurate  representation  of  the  actual  potential  for  a 
given  level  of  magnetic  activity.  Also,  if  it  can  be  de¬ 
termined  how  the  potentials  and  electric  fields  relate 
to  the  AE  indices,  then  it  will  be  possible  to  improve 
estimates  of  the  global  ionospheric  Joule  heating  as 
a  function  of  the  indices. 

1.1.  The  AE  indices 

The  Auroral  Electrojet  index  was  originally  intro¬ 
duced  by  Davis  and  Sugiura  (1966)  as  a  measure  of 
global  electrojet  activity  in  the  auroral  zone.  The 
index  is  derived  from  the  horizontal,  northern  com¬ 
ponent  of  geomagnetic  variations  observed  at  10  to 
13  observatories  along  the  auroral  zone  in  the  north¬ 
ern  hemisphere.  The  most  positive  value  measured 
from  all  the  stations  is  called  the  AU  index  and  the 
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most  negative  value  is  the  AL  index.  These  names 
refer  to  the  upper  and  lower  envelopes  of  the  over¬ 
lapped  plots  of  data  from  all  stations.  The  A I  and 
AL  indices  are  understood  to  provide  a  measure  of 
the  respective  eastward  and  westward  electrojet  cur¬ 
rents  at  a  given  time.  The  difference  between  the  AU 
and  AL  indices  gives  the  AE  index,  which  indicates 
the  total  maximum  amplitude  of  the  east  and  west 
electrojet.  currents.  The  mean  value  of  AU  and  AL 
is  the  AO  index. 

The  AU  index  is  usually  a  positive  number  and 
the  AL  index  is  usually  negative,  but  occasionally 
the  largest  horizontal  variation  (All)  is  a  negative 
number  and  at  other  times  the  smallest  variation 
(AL)  may  be  positive,  due  to  contamination  from 
the  magnetospheric  ring  current  or  zonal  currents  in 
the  ionosphere.  However,  the  AE  index  by  defini¬ 
tion  is  always  positive  and  is  independent  of  zonal 
currents,  although  the  physical  meaning  of  the  AE 
index  is  less  obvious  than  that  of  of  All  and  AL  in¬ 
dices  taken  separately  (Mayaud,  1980). 

The  term  “AE  indices”  is  commonly  used  to  re¬ 
fer  to  all  four  collective  indices.  These  indices  have 
advantages  over  the  standard  Kp  indices  in  that 
they  are  an  instantaneous  measurement  (i.e.,  one 
minute  rather  than  three  hour  resolution)  and  have 
interpretable  physical  meaning  (Davis  and  Sugiura, 
1966).  Presently  the  AE  indices  are  calculated  and 
distributed  by  the  World  Data  Center  C2  for  Geo¬ 
magnetism  at  Kyoto  University. 


1.2.  Throry 

At  this  point  it  is  prudent  to  consider  the  physics 
which  link  satellite  measurements  of  electric  poten¬ 
tial  to  the  AE  indices,  in  order  to  show  that  there 
is  a  valid  reason  for  such  a  comparison.  The  poten¬ 
tial  is  derived  by  an  integration  of  the  electric  field 
which  is  measured  by  the  satellite  along  it’s  direction 
of  motion: 


This  potential  is  relative  to  the  value  at  x  =  ib].  As¬ 
suming  that  the  measured  electric  field  is  meridional 
(north-south),  then  there  is  a  related  Hall  current  in 
the  azimuthal  (east-west)  direction: 


Jh  -  S/zEir  (2) 


where  is  the  height-integrated  Hall  conductivity 
and  J h  is  the  Hall  current  per  unit  length,  integrated 
over  the  thickness  of  the  ionosphere.  The  total  Hall 
current  is  given  by  the  integral  of  (2): 


hi r 

■fw. 


Ex  dr 


(3) 


where  E fj  has  been  assumed  to  be  constant  in  the 
integration  region.  If  u? j  and  ?/’•)  define  the  bound¬ 
aries  of  the  auroral  zone  (equatorward  of  the  convec¬ 
tion  reversal)  then  the  total  auroral  electrojet  cur¬ 
rent,  / // ,  is  proportional  to  the  potential  across  the 
auroral  zone: 

hi  £**  •  H) 

If  E//  is  not  constant  then  the  average  value  of  the 
conductance  in  the  region  may  be  substituted  as  an 
approximation. 

Since  the  AU  index  is  related  to  the  eastward  elec¬ 
trojet,  which  tends  to  be  located  in  the  evening  auro¬ 
ral  zone  (Davis  and  Sugiura,  19fi6)  and  is  a  Hall  cur¬ 
rent  driven  by  the  northward  electric  field,  then  the 
All  index  should  be  related  to  the  potential  across 
the  evening  aurora.  Likewise,  the  AL  index,  which 
corresponds  to  the  westward  electrojet,  should  be 
related  to  potential  across  the  morning  aurora. 
And  since  the  polar  cap  potential  is  equal  to  the  sum 
of  the  dusk  and  dawn  potentials,  and  the  AE  index 
is  equal  to  the  sum  of  AU  and  AL  (using  absolute 
values),  then  these  two  quantities  should  be  related. 
There  is  one  serious  problem  with  this  reasoning, 
however.  The  AU/AL  indices  are  not  a  measure¬ 
ment  of  the  total  electrojel  currents,  but  are  related 
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FIG  1.  SKETCH  OF  AN  ASSUMED  ELECTRIC  FIELD 
AND  ELECTROJET  CURRENT  DISTRIBUTION  IN  THE 
AURORAL  ZONES.  The  electric  field  has  a  triangular  distri¬ 
bution  with  peak  magnitude  E„  and  width  W.  The  Hall  (elec¬ 
trojet)  current  has  peak  magnitude  J(,  at  x  0.  The  electrojet 
current-  is  at,  altitude  above  the  ground,  where  it  causes  a 
horizontal  magnetic  perturbation  dBr.  The  magnitude  of  the 
greatest  horizontal  perturbation  is  obtained  by  integrating  the 
contribution  due  to  Jy(x)  in  the  interval  from  x  -W]  to  x:  w2. 
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lo  the  maximum  density  of  the  electrojet  current,  as 
measured  by  the  magnetic  observatory  which  hap¬ 
pens  to  be  closest  to  the  current  maxima  (Akasofu, 
1989). 

It  is  still  possible  to  derive  an  approximate  func¬ 
tional  relationship  between  the  potential  and  the 
peak  measured  ground  magnetic  perturbation  if 
some  assumptions  are  made  about  the  distribution  of 
the  electric  field  and  conductivity.  We  assume  that 
the  convection  electric  field  in  the  auroral  zone  can 
be  represented  by  a  triangular  function  with  width 
U  and  peak  E0,  as  shown  in  Figure  1,  and  the  Hall 
conductivity  is  constant.  The  integrated  electric  field 
is: 


<J> 


F.0\\ 

2 


(5) 


The  horizontal,  north-south  magnetic  field  pertur¬ 
bation  at  the  earth’s  surface  will  have  its  peak  value 
directly  below  the  location  of  the  maximum  Hall  cur¬ 
rent.  The  total  magnetic  perturbation  ran  be  calcu¬ 
lated  bv  integrating  the  contributions  from  infinite- 
hue  currents  of  magnitude  J y(r).  Using  the  well- 
known  formida  for  the  magnetic  field  around  an  in¬ 
finite  wire, 
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and  I  he  gpor^pfry  e}imvn  in  Figure  1,  the  total  hori¬ 
zontal  magnetic  perturbation  is 

Mo  fu’2  Jv(t)  cos9(x) 
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For  the  section  ot  the  electrojet  extending  from  x  -  0 
(the  location  of  the  peak)  to  x  =  w\,  the  integral  is: 


A  HT 


[L0J0h 
'2ir  j0 


rw i  1  -  -5- 

1  rfi* 


(8) 


where  ./„  is  the  peak  Hall  current  and  h  is  the  alti¬ 
tude  of  the  electro  jet  above  the  ground.  This  integral 
evaluates  to 

finJof(R\) 


a  nT 


2ir 


with  R i  and  /( R)  defined  as 

R I  7-  ^ 


(9) 


(10) 


f(R)  arctan  R 


2  R 


!n(  Rr  4  1).  (11) 


Adding  the  contribution  for  the  integral  from  0  to  w 2 
and  doubling  the  result  due  to  image  earth  currents, 
we  arrive  at 

1  f{R  2)1 
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Using  (2)  and  (5).  the  final  result  is 


11  7r 

Since  the  altitude  of  the  electrojet,  h,  is  about  100 
km,  and  the  half-widths  are  in  the  range  of  250  to 
1,000  km,  then  R  is  in  the  range  of  2.5  to  10.  and 
f(R)  varies  from  0.79  to  1.2-1;  the  function  /  is  not 
very  sensitive  to  the  width  or  altitude  of  the  electro- 
jet.  Equation  13  ran  be  reduced  to  the  form 

AHX  ■-  c£/y4>  .  (14) 

I  he  factor  c  is  expected  to  vary  from  approximately 
0.99  nT-ohtn-kV  *  for  W— 2000  km  to  2.54  nT-ohtn- 
kV  1  for  W  -  500  km. 

In  this  analysis  the  Hall  currents  in  the  polar  rap. 
which  would  reduce  the  magnitude  of  A BT.  have 
been  neglected  as  these  currents  are  usually  too  far 
away  to  have  a  significant  effect.  However,  there  is 
one  more  aspect  to  this  problem  which  has  more  sc 
rious  consequences:  the  electric  field  measured  by 
the  satellite  (Equation  1)  is  in  a  co-rotating  refer¬ 
ence  frame,  while  the  electric  field  which  determines 
the  ionospheric  currents  (Equation  2)  is  in  a  frame 
relative  to  the  neutral  wind  velocity.  The  neutrals 
are  usually  accelerated  in  the  direction  of  E  x  B 
by  drag  forces  (Killeen  et.  al..  1984).  This  has  the 
effect  of  making  the  electrojet  current  smaller  than 
what  would  be  calculated  from  the  measured  poten¬ 
tials  and  the  expected  Hall  conductivity.  Lyons  and 
Waltersheid  (1986)  have  estimated  that  the  neutrals 
will  attain  speeds  of  approximately  40%  of  the  ion 
convection  velocity,  which  would  result  in  an  “effec¬ 
tive  flail  conductivity”  which  is  60%  of  the  actual 
Hall  conductivity.  Simultaneous  measurements  of 
ion-drift  velocities  and  neutral  wind  velocities  show 
that  the  relation  between  the  two  is  complex  and 
highly  variable.  Killeen  r(  al.  (198-1)  have  found 
that,  at  times  the  neutral  and  ion  velocity  vectors 
show  agreement  in  both  magnitude  and  direction, 
and  at  other  times  or  locations  there  can  be  signif¬ 
icant  differences.  They  also  found  that  the  veloci¬ 
ties  in  the  dusk  auroral  zone  tended  to  match  better 
than  in  the  dawn  convection  channel.  The  result  of 
this  complex  ion-neutral  relationship  is  to  introduce 
additional  variability  into  the  relation  between  the 
magnetic  field  perturbations  and  the  integrated  po¬ 
tentials. 

We  also  need  to  keep  in  mind  that  the  actual  mea¬ 
surements  are  subject  to  some  errors  which  may  up¬ 
set  the  expected  correlation.  When  the  auroral  oval 
contracts  poleward  at  times  of  low  activity  many  of 
the  AE  observatories  will  be  well  equatorward  of  the 
electrojet  currents  and  will  not  make  reliable  mca 
surements.  And  at  times  the  satellite  may  measure 
a  fraction  of  the  full  potential,  since  the  orbit  path 
might  not  cross  the  regions  where  the  electric  poten- 
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tial  reaches  its  minimum  and  maximum  values. 

1.3.  Previous  publications 

The  theoretical  relationship  between  the  polar  cap 
potential  and  ground  magnetic  perturbations  had 
been  considered  by  Crooker  and  Siscoe  (1981),  who 
used  a  two-ring  model  of  the  auroral  zone  and  polar 
cap.  However,  their  expression  for  the  magnetic  per¬ 
turbations  is  given  for  a  position  equatorward  of  the 
auroral  oval,  rather  than  the  maximum  values  (AE) 
within  the  oval.  Right  at  the  equatorward  boundary, 
their  Equation  A 17  reduces  tor 

±B  *  (15) 

■lb* 

where  4>  is  the  total  polar  cap  potential,  a  is  the 
radius  of  the  polar  cap  (inner  radius  of  the  auroral 
oval)  and  h  is  the  outer  radius  of  the  auroral  zone. 

Another  theoretical  analysis  has  been  done  by  Nis- 
bet  (1982),  in  which  the  relations  between  the  Birke- 
land  currents.  Joule  heating,  and  the  AE  indices  are 
examined.  Nisbet’s  Equation  3d  for  the  relationship 
between  the  electric  potential  and  A  Pi  indices  is  iden¬ 
tical  in  form  to  our  Equation  13.  However,  his  ex¬ 
pression  (Equation  23)  for  the  factor  /  is  different, 
even  though  it  is  also  based  on  a  triangular  distri 
bution.  (It  appears  that  the  cos  8  adjustment  was 
omitted.)  Nisbet  did  note  that,  since  the  polar  cap 
potential  shows  less  variation  than  the  AL  and  AU 
indices,  the  variations  in  the  electrojet  indices  are 
“mainly  a  response  to  Hall  conductivity  variations 
rather  than  to  variations  in  the  generator  driving  the 
Birkeland  currents”. 

On  the  experimental  side,  the  publications  related 
to  the  present  topic  are  much  more  numerous.  Satel¬ 
lite  measurements  of  elect  ric  potential  across  the  po¬ 
lar  rap  are  compared  to  solar  wind  parameters  by 
ReifT  et.  al.  (1981),  Wvgant  el  al.  (1983),  and  Doyle 
and  Burke  (1983).  ReifT  ci  al.  (1981)  had  found  that 
among  the  various  indices  of  magnetic  activity,  'tie 
AE  index  had  the  highest  correlation  with  polar  can 
potential,  although  the  linear  correlation  was  not  as 
good  as  with  the  solar  wind  parameters.  The  best 
lit  was  found  to  be  <t>(kV)  41  +  0.1  lAE(nT).  Since 

their  primary  purpose  was  to  investigate  the  solar- 
wind  coupling,  this  relation  was  not  explored  in  any 
detail. 

There  have  also  been  studies  of  the  relationship 
between  solar  wind  parameters  and  the  AE  indices 
by  Rerrault  and  Akasofu  (1978),  Meng  et  al.  (1973), 
Baker  et  al.  (i981),  Bargatze  et  al.  (1985),  Tsuru 
tani  el  al.  (1985),  and  Sauvaud  et  al.  (1987)  (to 
name  just  a  few).  A  common  element  in  several  of 
these  papers  is  'he  finding  that  the  best  correlation 
between  the  AL  index  and  IMF  is  found  when  the 
AL  index  is  lagged  by  40-60  minutes. 


Various  models  for  the  distribution  of  the  electric 
potential  or  plasma  convection  over  the  polar  cap 
have  been  formulated,  but  these  models  have  not 
been  related  to  the  AE  indices.  Empirical  convection 
patterns  have  been  derived  by  Foster  et  al.  (1986) 
from  radar  observations:  these  patterns  are  keyed  to 
"precipitation  indices”  which  were  derived  from  par¬ 
ticle  measurements  on  NOAA/TIROS  weat  her  sat  el 
lites.  A  linear  relationship  was  found  between  the  po¬ 
lar  cap  potential  and  the  precipitation  index,  which 
varies  from  1  to  9.  In  another  publication,  Foster  ct 
ill.  (1989)  have  also  calculated  ionospheric  conduc¬ 
tances  and  field  aligned  currents  for  the  same  precip¬ 
itation  indices. 

Ileppner  and  Maynard  (1987)  had  derived  empiri¬ 
cal  convection  models  from  DE-2  electric  field  flat  a. 
keyed  to  various  signs  of  z  and  ;/  components  of  the 
interplanetary  magnetic  field.  These  models  were 
drawn  for  moderately  active  conditions  (Kp  3'  ). 
Average  cross  polar  cap  potentials  for  different  lev¬ 
els  of  Kp  and  AE  were  given  and  the  expansion  of 
the  lower  boundary  with  Kp  was  shown.  Rich  and 
Maynard  (1989)  put  the  Ileppner -Maynard  models 
into  analytic  form,  with  the  magnitude  of  the  po¬ 
lar  cap  potential  depending  oil  the  Kp  index.  Iono¬ 
spheric  currents,  joule  heating,  and  field  aligned  cur 
rents  were  then  calculated,  using  the  analytic  electric 
field  models  in  combination  with  ionospheric  conduc¬ 
tivities  which  were  based  on  a  precipitation  model  by 
Hardy  rt  al.  (1987).  This  precipitation  model  is  also 
based  on  Kp. 

Spuo  "t  al.  (1982)  have  derived  models  for  precipi¬ 
tating  electron  energy  and  auroral  zone  conductances 
which  are  parameterized  according  to  the  AE  index, 
it  was  found  that  the  AE  index  does  a  better  job  of 
ordering  the  data  than  other  indices  of  geomagnetic 
activity.  It,  would  be  useful  to  know  how  the  elect  ric 
potential  varies  according  to  the  same  index. 

There  have  been  many  publications  by  Y.  Kamide 
and  associates  about  a  method  for  calculating  the 
distribution  of  currents  and  electric  fields  over  the 
entire  polar  rap  and  auroral  oval,  through  an  “in 
version”  of  the  magnetic  records  from  a  number  of 
stations  in  the  northern  hemisphere.  Kamide  and 
Baumjohann  (1985)  show  both  the  total  estimated 
electrojet  currents  and  potential  differences  com¬ 
pared  to  the  AU/AL  indices.  The  ratio  between 
the  total  westward  electrojet.  and  AL  shows  more 
variation  than  the  eastward  electro jet  AT  ratio,  par 
tirularly  during  substorms.  They  also  observe  that 
the  ratio  of  the  total  electrojet  current  to  the  to¬ 
tal  potential  difference  increases  during  substorms. 
Kamide  and  Baumjohann  (1985)  conclude  that  the 
increase  in  the  ionospheric  current  intensity  prior  to 
substorm  onsets  is  due  to  an  increase  in  the  electric 
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fit-id.  and  conductivity  enhancements  aftt’r  onsets 
contribute  to  the  maximum  sub.itorm  currents.  Alin 
ft  id.  (1989)  make  the  same  conclusion,  also  based 
on  the  Kamide  method  of  inserting  the  ground-based 
magnetograms. 

Ahn  rf  nl.  (1981)  have  also  used  the  inversion 
technique  to  compute  polar  cap  potentials  contin¬ 
uously  for  a  two-dav  period  in  March,  1978.  The 
derived  potentials  are  then  compared  to  the  simulta¬ 
neous  AE  indices.  Their  linear  regression  resulted 
in  <J>(kY)  -  36  +  0.089AE] 2(nT)  when  12-station 

rnagnetograms  were  used  in  the  inversion.  With 
more  accurate  71-station  computations  the  result 
was  4>(kV)  2-1  *  0.098AEy](nT).  Baumjohann  and 

Kamide  (1984)  used  the  same  technique  to  derive 
•i  liiivu.  i egression  between  total  ionospheric  Joule 
heating  and  the  AF  indices.  However,  their  data 
shows  a  definite  non-linear  trend,  which  they  did  not 
fit  due  to  the  sparsity  of  data  for  high  AE  values. 
More  recently,  Richmond  ci  al.  (1990)  have  used  hr 
inversion  technique  combined  with  incoherent  scatter 
radar  observations  from  a  two-day  period  in  January, 
1981.  They  found  that  4>(kV)  -  22 -t  0.1 19AEj2(nT). 

Ahn  cl  al.  (1983)  derived  empirical  relations  be¬ 
tween  ground  magnetic  disturbances  and  local  iono¬ 
spheric  conductivities  and  electric,  fields,  based  on  a 
comparison  of  data  from  one  magnetic  observatory 
and  a  nearby  radar  facility.  The  electric  fields  were 
related  to  the  magnetic  perturbations  by  power  law- 
equations,  with  different  coefficients  for  positive  and 
negative  magnetic  disturbances. 

We've  seen  various  parameters  which  have  been 
related  to  the  AE  indices  for  limited  time  intervals. 
However,  it  is  well  known  that  the  ionospheric  con¬ 
ductivity  changes  with  the  season.  Burch  (1973) 
found  seasonal  variations  in  AU,  but  not  in  AL. 
Burch  noted  that  this  was  consistent  with  the  results 
of  the  conjugate  point  study  of  Meng  and  Akasofu 
(1908)  which  showed  that  the  positive  bays  in  the 
summer  hemisphere  were  roughly  twice  those  in  the 
winter  hemisphere,  but  there  was  no  clear  asymme¬ 
try  for  negative  bays.  Allen  and  Kroehl  (1975)  also 
found  significant,  systematic  seasonal  AU  intensity 
differences  which  average  60  nT  between  winter  and 
summer.  The  AL  values  were  not  as  clearly  differen¬ 
tiated.  Berthelier  (1976)  also  saw  seasonal  effects  in 
A  F,  which  were  attributed  to  conductivity  variations. 
There  were  also  annual  and  diurnal  variations  in  AE 
which  Berthelier  attributed  to  a  shift  of  B^(GSM) 
relative  to  B^(CJSE),  which  has  a  symmetrical  dis¬ 
tribution  around  a  zero  mean. 

This  synopsis  of  the  literature  indicates  that  a 
comprehensive,  experimental  comparison  of  AE  in¬ 
dices  and  direct  measurements  of  potential  is  still 
larking.  Furthermore,  seasonal  variations  in  iono¬ 


spheric  conductivity  nerd  to  lit-  taken  inlo  consider 
a  l  ion. 

2.  OBSERVATIONS 

2.1.  J'hr  DK-2  data  set 

The  DE-2  satellite  was  launched  in  August  of  1981 
into  a  300  to  1009  km  altitude  polar  orbit.  The  satel¬ 
lite  ree:  ’ered  the  atmosphere  in  March  of  1983.  Elec¬ 
tric  fields  were  measured  on  DE-2  by  two  orthogo¬ 
nal  double  probes  in  the  orbit  plane,  with  lengths 
of  22.1  m  tip-to-tip.  The  data  rate  was  16  samples 
per  second.  Further  details  about  the  Vector  Elec¬ 
tric  Field  Instrumen  (VEFI)  are  given  by  Maynard 
rt  al  (1981). 

Electric  potentials  are  obtained  by  an  integra 
lion  of  the  measured  electric  fields  along  the  satel 
lite’s  path.  For  this  study  VEFI  data  were  analyzed 
from  309  auroral /polar  cap  passes  in  the  period  from 
November  6,  1981  to  January  4,  1982;  and  218  passes 
in  the  period  from  May  1,  1982  to  June  29.  1982.  The 
data  in  the  first  period  correspond  to  ..inter  con¬ 
ditions  in  the  northern  hemisphere,  and  the  second 
period  correspond  to  summer  conditions.  The  mag¬ 
netic  local  times  of  the  satellite’s  orbital  track  during 
these  time  periods  was  within  about  3  hours  of  dawn 
and  dusk  (6  and  18  hours  MLT).  This  restriction 
is  given  so  that  the  satellite  will  have  a  reasonable 
chance  of  passing  through  the  potential  maximum 
and  minimum  of  the  typical  two-cell  convection  pat¬ 
tern.  The  data  used  in  this  study  do  not  include  all 
that  was  obtained  by  the  VEFI  instrument  during 
these  t;  ne  periods.  There  was  some  random  selec¬ 
tion  of  events  due  to  circumstances  of  availability  of 
processed  data,  and  some  passes  were  unsuitable  for 
a  measurement  of  the  polar  cap  potential. 

Electric  fields  measured  with  double  probes  are  not 
without  some  uncertainty.  “Contact  potential”  dif¬ 
ferences  on  the  two  probes  can  cause  a  small  offset,  in 
the  measured  electric  fields,  which  can  lead  to  large 
errors  in  the  potentials  when  integrated  over  long  dis¬ 
tances.  To  minimized  these  errors  the  contact  poten¬ 
tials  were  adjusted  for  every  polar  pass  in  this  study. 
Adjustments  were  made  in  such  a  way  that  the  dif¬ 
ference  in  integrated  potential  between  the  points  of 
the  pass  at  45  degrees  invariant  latitude  on  each  side 
of  the  pole  was  zero.  In  some  cases  the  pass  (when 
the  instrument  was  turned  on)  started  or  ended  at 
latitudes  higher  than  45  degrees.  In  those  cases  the 
beginning  or  end  of  the  pass  was  used  for  the  adjust¬ 
ment,  provided  that  the  electric  field  was  still  flat  at 
that  invariant  latitude.  If  the  beginning  or  the  end 
of  the  pass  was  already  inside  the  auroral  zone  <  ,‘n 
the  pass  was  not  used  for  this  study. 

The  high  latitude  boundaries  of  the  dawn  and  dusk 
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(a)  (b) 

Kir,.  2.  LOCATION  OF  THE  PASSES  BY  THE  DR-2  SATELLITE,  IN  MAGNETIC  LOCAL  TIME  • 
INVARIANT  LATITUDE  COORDINATES,  FOR  THE  PERIOD  FROM  6  NOVEMBER.  1981  THROUGH 
4  JANUARY,  1982.  (a)  Location  of  the  dawn  and  dusk  auroral  zone  passes,  (bj  Location  of  the  polar  cap 
passes 


FIG.  3.  LOCATION  OF  THE  PASSES  BY  THE  DE-2  SATELLITE,  IN  MAGNETIC  LOCAL  TIME  - 
INVARIANT  LATITUDE  COORDINATES,  FOR  THE  PERIOD  FROM  1  MAY,  1982  THROUGH  29  JUNE. 
1982.  (a)  Location  of  the  dawn  and  dusk  auroral  zone  passes,  (b)  Location  of  the  p",r*r  passes. 


I'okii  cap  potentials  ami  auroral  electrojet  indices 


auroral  zones  were  determined  by  the  local  maxima 
and  minima  of  the  integrated  potential  curves.  The 
!ct  latitude  boundaries  of  the  auroral  zones  were  de¬ 
termined  bv  choosing  the  points  in  which  the  electric 
field  curve  stopped  being  flat  and  rapid  variations 
started.  Figure  2a  shows  the  locations  of  the  dawn 
and  dusk  auroral  zones  derived  from  the  electric  field 
data  for  the  winter  period.  This  graph  shows  the 
orbit  tracks  of  the  DE-2  satellite  between  the  low 
and  high  latitude  auroral  boundaries,  in  invariant 
latitude-MLT  coordinates.  In  Figure  2b  are  shown 
the  polar  cap  portions  of  the  same  passes.  Figures  3a 
and  3b  have  the  same  format,  but  show  the  locations 
of  the  measurements  during  the  summer  period. 

The  difference  between  the  potentials  at  the  local 
maxima  and  minima  gives  the  total  polar  cap  po¬ 
tential.  The  difference  between  the  potentials  at  the 
low  and  high  latitude  boundaries  of  the  auroral  zones 
gives  the  total  potential  across  the  dawn  and  dusk 
auroral  zones.  These  relationships  are  illustrated  it. 
Figure  4,  which  shows  an  example  of  the  electric  field 
over  the  northern  polar  cap  and  the  associated  po¬ 
tentials. 

2.2.  Comparison  with  AE  Indices 

The  AE  indices  with  one  minute  time  resolution 
were  obtained  on  magnetic  tape  directly  from  the 
World  Data  Center  C2  for  Geomagnetism  at  Kyoto 
University.  An  average  AU  index  was  obtained  for 
each  dusk  auroral  zone  crossing  by  averaging  the  one- 
minute  AU  indices  corresponding  to  just  the  time 
period  that  the  DE-2  satellite  was  in  the  dusk  au¬ 
roral  zone,  as  determined  by  the  electric  field  sig¬ 
natures.  Likewise,  an  AL  index  was  determined  for 
each  dawn  auroral  zone  crossing,  and  an  AE  index 
was  determined  for  each  polar  cap  crossing.  Thus, 
our  comparisons  are  between  simultaneous  measure¬ 
ments. 

Scatter  plots  of  AU  indices  which  are  plotted  in 
comparison  with  the  dusk  auroral  potentials  are 
shown  in  Figure  5.  The  signs  of  the  dusk  poten¬ 
tials,  which  are  shown  as  negative  numbers  in  Fig¬ 
ure  4,  have  been  reversed  on  this  graph.  There  are 
separate  plots  for  each  of  the  two  time  periods.  As 
mentioned  earlier,  one  time  period  corresponds  to 
winter  in  northern  hemisphere  and  the  other  time 
period  corresponds  to  summer  in  the  northern  hemi¬ 
sphere.  The  measurements  which  were  made  while 
DE-2  was  in  the  northern  hemisphere  are  marked 
with  the  crosses,  and  southern  hemisphere  passes  are 
marked  with  the  diamonds.  Therefore  the  symbolic 
representation  of  the  summer/ winter  “season”  is  re¬ 
versed  beiween  the  two  plots,  but  the  crosses  consis¬ 
tently  show  potential  measurements  in  the  northern 
hemisphere,  where  the  ground  magnetic  observato¬ 
ries  are  located. 


There  are  three  dashed  lines  drawn  on  each  of 
these  graphs.  The  lowermost  dashed  line  shows  a 
linear  relationship  between  AU  and  the  dusk  poten¬ 
tial,  with  the  slope  chosen  such  that  10%  of  the  data 
points  lie  below  the  line  and  90%  of  the  points  are 
above.  Likewise,  the  uppermost  dashed  line  shows 
the  slope  at  which  10%  of  the  points  are  above  the 
line  and  90%  are  below  These  slopes  bound  the  most 
likely  range  of  the  product  “cE^”  in  Equation  14. 
The  middle  dashed  line  shows  the  median  slope,  at 
which  50%  of  the  points  are  above  the  line  and  50% 
are  below  the  line.  The  values  of  the  10%,  50%,  and 
90%  slopes  are  given  in  Table  1. 

In  Figure  6  there  are  shown  similar  scatter  plots  for 
the  AL  indices  and  the  simultaneous  measurements 
of  the  dawn  auroral  potential.  The  signs  of  the  AL 
indices  have  been  reversed.  Lines  showing  I  he  10%, 
50%,  and  90%  slopes  are  also  indicated  on  the  figure. 
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FIG.  4.  EXAMPLE  OF  ELECTRIC  FIELDS  MEASURED 
BY  THE  VEFI  INSTRUMENT  ON  DE-2  (upper  panel)  AND 
POTENTIALS  OBTAINED  BY  INTEGRATING  THE  ELEC¬ 
TRIC  FIELD  (lower  panel).  The  dusk  and  dawn  potentials  are 
defined  as  the  difference  between  the  potentials  at  the  low  and 
high  latitude  boundaries  of  the  respective  dusk  and  dawn  au¬ 
roral  sones.  The  polar  cap  potential  is  defined  as  the  difference 
between  the  potentials  at  the  two  high  latitude  boundaries, 
which  are  determined  by  the  local  minimum  and  maximum. 
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FIG.  5.  SCATTER  PLOT  OF  AU  INDEX  VS.  DUSK  POTENTIAL,  (a)  Winter  period  in  the  northern 
hemisphere,  (b)  Summer  period.  The  crosses  mark  the  cases  where  the  potential  was  measured  in  the  northern 
hemisphere,  and  the  diamonds  show  the  southern  hemisphere  measurements.  The  dashed  lines  are  at  constant 
slopes  which  bound  10%,  50%.  and  90%  of  the  data  points.  The  parameters  for  these  lines  are  given  in  Table  1 
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FIG.  6.  SCATTER  PLOT  OF  AL  INDEX  VS.  DAWN  POTENTIAL,  (a)  Winter  period  in  the  northern 
hemisphere,  (b)  Summer  period.  The  format  is  the  same  as  in  Figure  5. 


TABLE  1  SLOPES  OF  THE  DASHED  LINES  ON  FIGURES  5,  6,  AND  7 


Data 

Season 

Figure 

10%  Slope 

50%  Slope 

90%  Slope 

AU  vs.  4>Du9k 

winter 

5a 

0.89 

2.07 

4.75 

summer 

5b 

1.61 

3.86 

8.38 

AL  vs.  $Dawn 

winter 

6a 

0.43 

3.65 

13.13 

summer 

6b 

2.88 

7.63 

4 1 .84 

AE  vs.  $[1(2 

winter 

7a 

1.17 

2.81 

6.52 

summer 

7b 

9  9/1 

4.95 

9.62 
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(a)  (b) 

FI<J.  7  SCATTER  PLOT  OF  AE  INDEX  VS  I’OI.AR  CAP  POTENTIAL  in  Winter  period  in  the  northern 
hemisphere  f  hi  Summer  period.  The  format  is  tlie  same  .us  in  Figure  r> 


(a)  (»>) 

FKI  8.  SCAT  TER  PLOT  OF  AE  INDEX  VS  POLAR  CAP  POTENTIAL.  WI  TH  THE  MEASUREMENT 
OF  AE  LAGGED  50  MINUTES  BEHIND  THE  MIDDLE  OF  THE  POLAR  CAP  CROSSING  (a)  Winter 
period  in  the  northern  hemisphere,  (b)  Summer  period. 


and  the  values  arc  included  in  Tabic  1 .  There  arc  sev¬ 
eral  points  in  Figure  6  which  lie  substantially  above 
and  to  the  left  of  the  90%  slopes.  Most  likely  this 
is  because  the  satellite  only  measured  a  fraction  of 
the  full  potential.  Most  of  these  abnormal  potentials 
were  measured  in  the  winter  hemisphere  (plus  sym¬ 
bol  in  Figure  da  and  diamond  symbols  in  Figure  Sb), 
and  at  magnetic  local  times  of  8  to  9  hours.  At  other 
local  times  there  are  no  significant  differences  be¬ 
tween  the  potentials  in  the  northern  and  southern 
hemisphejes. 


Figure  7  shows  the  results  of  comparing  the  A  F 
indices  to  the  polar  cap  potentials,  with  the  same 
format  as  in  F  igures  5  and  6.  Similarly,  (lie  constants 
which  describe  the  superimposed  lines  are  in  Fable  I . 
The  change  in  the  ordinate  axis  is  due  to  the  fact  I  hat 
AF,  is  by  definition  always  positive,  whereas  AF  and 
\L  may  have  either  sign. 

We  emphasize  that  the  pola.  rap  potentials  and 
AE  indices  which  arc  shown  in  Figure  7  are  sirnul 
taneous.  It  had  previously  been  shown  that  the  po¬ 
lar  cap  potential  is  related  to  the  IMF’  (Rciff  ct  til, 
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f't^j  i.  and  ••  had  also  been  shown  that  t  he  best  <  or 
relation  between  I  fit'  IMI*  and  A  I  is  obtained  when 
\  I  h-  lagged  by  M>  i«»  hit  minute-  (  Maker  »  t  >ii.  f  ‘IK i  ; 
Mnrga!  7v  */  ill  .  T  ?  I K  r»  I  I  herefore,  w-  have  asked  the 
r j i j < ■  s t  i i . i *  about  wIh'U'jt  or  n*t  a  delav  of  .'At I  min 
ii* » '=-  b«  tween  the  [>o!;ir  cap  potential  measurements 
and  t|n-  \  V  iru'iisiirrmcnis  would  reduce  t !i<*  scatter 
in  » f m ■  data  points.  I  hr  answer,  shown  in  Kigur«* 

•  s  ’hat  the  -ratter  is  made  worst’  by  lagging  the  \ !' 
measurement  behind  the  potential  measurement  by 
AM  rnonitf1' 

>in(*'  geomagnetic  Mjb-forms  by  definition  do 
u!M'  l,»  r  gr  <b-  \  ia  i  ioiii>  in  the  \  I .  m<U*x .  tin*  next  st  ep 
i-  to  determine  if  the  location  of  data  points  on  the 
-rutin*  graphs  is  related  to  substorm  phase.  Vor  all 
A'27  ra-rs  us^d  in  this  -tudv  the  charts  of  the  A  I'. 

huh  air  published  by  the  World  I  bit  a  <  enter  ( '2 
for  ( i  cor  nagne  t  i  Mu.  were  examined  to  see  if  the  satel¬ 
lite  measurements  1  ,M*k  place  during  a  well  defined 
s'ibs»..nii  -growl h  phase",  “expansion  phase".  >*r  “re 
co’.er\  phase".  In  main  cases  an  unambiguous  phase 
n>uhi  riot  be  identified.  for  *  he  points  where  the 
phase  was  <  lear.  the  results  are  shown  in  Vigor***  M 
and  H)  I  he  results  indicate'  that  the  substorm  phase 
doe**  indeed  influence  the  location  of  pmnts  on  the 
graphs 

3.  discission 

I  h«-  re- ui'  '  simw  n  heir  nidi*  a*e  *  hat.  as  evpec t*  r|. 

!  f><  rr  la  1 1*  «n  **h  i  ps  between  the  \  V  indues  and  t  fie 

*  n*  a s i j re <  1  auroral  polar  cap  potentials  depend*  «-n 
ih<  time  *>f  the  vc;i i .  I  he  A  I*  indices  are  larger  in 
magnitude  when  it  is  summer  in  the  northern  ii'*mi 

-  pfiere .  where  thev  are  measured.  due  to  the  greater 
i"nospherit  (  omf  ip  t  i  vii  v  i  Allen  and  Kroehl.  IM7Ai 
I  hi"  i*  evident  in  ttie  upward  shift  m  the  distribution 
of  the  \  f-  indices  between  1  igures  7a  and  7b  I  here 
was  abo  a  downward  shift  in  the  distribution  of  the 
measured  polar  cap  potential  between  t  lie  two  lime 
periods,  which  can  be  attributed  to  a  lowei  sunspot 
number  in  June  10^2  compared  to  December  l!J*l. 
f  or  the  rases  where  IMI*  measurements  arc  available 
we  have  also  fori  nr  j  that  the  di  st  ri  bn  t  i<»n  of  It  ( J  S  M  ) 
had  *  hanged  In  more  northward,  which  accounts  for 
t  lie  decreased  polar  cap  potentials. 

In  the  graphs  which  have  been  shown  here  the  data 
have  been  separated  according  to  the  season  in  the 
northern  hemisphere,  where  the  magnetic  observalo 
rif*s  are  located,  fins  give^  a  mom  uniform  group 
ipg  of  flu*  data  points,  in  comparison  to  graphs  (not 
sh*  >w  u  )  w  hit  h  group  the  tbit  a  ac»  ording  to  t  he  seas*  mi 
local  to  the  satellite  "  potential  measurement*  li  e., 

combining  the  November  January,  northern  lemi 
sphere  passes  with  the  May  June,  southern  hemi 
sphere  passe*  |. 


It  i*  expected  that  til**  electric  potential  -ho'iid 
be  the  same  over  the  northern  and  southern  polar 
laps,  since  thep  boundaries  are  connected  b;  t  he 
<  lunluet  ing  magneti*  field  lines.  I  he  data  do  pot 
disagree  with  this  assumption,  -mr**  tin-  southern 
hemisphere  measurement-  tend  to  relate  to  tin  \l 
indices  in  the  same  way  as  the  northern  hepuspher* 
measurements.  Ib*wever.  there  ar*  cases  v*.  her*  t  h* 
pi  m  e  n  *  in !  iiie,iMir  emerii  *  over  the  winter  h»  ttu«-phe» i- 
seem  t*>  be  less  accurate  than  the  summer  mea-m* 
merit  s  It  a  p  pears  that  the  f  i  \  n*  h  n  a !  i*  if:*  "f  1  }  i#-  p<  *  j  o  : 
cap  auroral  boundaries  tend  to  be  ohw  ured  bv  *  ** 
erratic  elect  ric  field  fluci  ital  ions  w  hh  h  are  as  so*  m1  *  ■; 
with  h«w  i< urns phene  eond net  i vi » y .  pa r !  lcniari  \  i r  ‘  - 
late  morning  auroral  /on r.  I  lie  iru  r*  a-»*ri  ero  * 
the  winter  is  fi  om  1  he  :na  bi  I  i  tv  f  <  ■  -  a  - 1 1  •.  !*!*■•■•'' 
reversal  boundary. 

I  lie  ele<  I  ro  |et  *  line  lit  s  an*  I  t  hei  r  *  ‘  •  *  *•-  ,o  <  >ij. 
to  be  approximately  pt  <*pr -r  t  i<-na!  to  tn*_  pof<u;''.« 
as  predicted  by  l.quaiiofi  I  I  I  fir**  i-  ,i  fair  ;«p!>>  ■  ’ 
of  spread  !<>  the  rlata  **ifiM  tin-  wid;h  and  •  i  •  -  *  *  1  > : . 

J  ion  of  the  r  !*•*  I  r*  >  at  <  ;  j  ri  *  r  *  and  t  li«  1 1  a !  1  <  «  up!  ■  : 

1 1  v  are  not  constant,  f  *  •  r  a  i«unparisori  b*lw**":  ■  :• 
data  and  .  |h..h,  hi  u-  u-e  v.iluc*.  of  solar  rad*a 1  ;■ 
induc«*d  «ondu*ti\ii\  erpial  *  *  *  n  in  tin  wirite*  .-vi 
7  mhos  in  the  summer  \Ne  need  ?*•  and  t<>  tins  b,p  k 
gr**und  ( ondiict  i  vity  i  taking  the  sp.,.i?e  r<*ot  <d  •  h* 

s  1 1  f  r  1  of  tile  - 1 J 1 1  a  1*  -S  1  I .  *  •  <  o  n  *  1 1 !  *  *  I  ■  I  *  ■  I  *  w  h  M  *  •  *  I  1  >  •.  n  1  ' 

fi*le  pre<  ipit  at  ion.  \i  cording  ffard\  ■(  d 
f»  *  r  K.,  '.}  tins  induced  c**nduct  \  vity  i<  a  pproxitM.t  •  r!  \ 

J  mho  at  1  s  hours  M  [  I  ami  !  t»  mhos  at  b  fioiirs  M  !  I 
I  herefore  in  the  dusk  electrojet  the  r. *mluct i vit y  i*l 
approximate! v  J  mh*»  in  winter  and  7  fi  mfe*  in  s*m: 

I  per.  an*  I  the  slope's  on  \  igtire  A  sfumld  be  about  d  7"t 
and  It  A  n»  k\  ;  In  t  he  da  wn  rlert  r<*  jet  the  condur 
livitv  is  about  III  mho  m  winter  and  1  ‘1  '1  tnh»*  !*■■ 
summer,  and  tin*  slopes  *mi  l  igure  b  should  b*>  a'>out 
12  A  and  I  A. 2  til  k\  1  We  find  that  mir  ‘dM'A  **:»>]'.-- 
whicli  are  give*n  in  I  able  I  ;ir<*  not  far  from  these 
values,  but  <ui  the  average  the  slopes  are  much  !e-s 
Most  likely  this  is  due  f<>  the  effects  <»f  the  neutral 
winds,  as  discussed  previously. 

I  tip*  to  these  numerous  other  variables  which  v\e 
cannot  control,  the  data  show  a  large  amount  of 
'scatter  on  the  graphs.  Nevertheless,  there  is  suffi 
*  h  id  older  to  enable  us  to  fit  t  he  data  t«*  a  curve  An 
equation  which  describes  the  relationship  between 
t  fie  A  I*  index  and  polar  cap  potential  drop  would 
have  some  utility,  as  mentioned  in  the  In*  roduct  ion. 

\\e  will  concentrate  on  just  the  data  shown  in  Vig 
nre  7,  but  first  we  will  attempt  t:.  eliminate  from 
the  data  base  some  of  the  satellite  passes  which  did 
not  pass  through  the  full  potential  drop.  H  is  well 
known  that  the  plasma  convection  pattern  tends  to 
be  rotated  slight  I  v  clockwise  (  Hep  p  tier  and  Maynard. 


} !  i ;  1 

I  II!-’  ‘ 
I  HI  t 
I’M  \S 


I't'i.ii  and  .iunu.ll  cka  i:  i  •  1 1  !  until  - 


/  AE  INDEX  \S  I't  >1  AH  CAP  POTENTIAL  FOR 
SOUTHERN  H  KM  IS  I’M  FR  E  WIN  ITR  PERIOD.  FOR 
ASI's  IN  WHICH  AN  I  N  A  M  RK  HOPS  SP  BSTORM 
I  CUi'l.ll  RF  IDENTIFIED  Hirst*  data  points  were 
r*H  ri'jrinn  '  ^rn’Aih,  ih‘  expansion.  and  <  r ! 

’  •.  js!  ;n»f  s 
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rrnHciv  stages. 
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1987),  and  in  Figure  6  it  was  found  that  some  of 
t fie  dawn  potentials  which  were  measured  at  MLT 
greater  than  8  hours  were  obviously  low.  There¬ 
fore,  we  use  just  the  orbits  wh:ch  passed  through  the 
dawn  aurora!  zont  at  MLT  less  than  8  hours,  and 
also  passed  through  the  dusk  auroral  zone  at  MLT 
less  than  20  hours.  This  reduces  the  number  of  mea¬ 
surements  from  309  to  211  cases  in  the  winter  time 
period,  and  from  218  to  129  cases  in  the  summer 
time  period. 

Conventional  least-squa .e-ertor  techniques  could 
he  used  to  determine  the  curve  fits.  However,  these 
methods  minimize  the  vertical  error  between  the 
curve  and  the  “dependent  variable”  data  points,  as¬ 
suming  that  the  “independent  variable”  on  the  hor¬ 
izontal  axis  has  been  measured  with  a  much  smaller 
error  A  least  square-error  fit  with  AE  as  the  inde 
pendent  variable  gives  a  result  which  is  substantially 
ilillerent  from  the  result  which  is  obtained  with  $, 
the  polar  cap  potential,  as  the  independent  variable. 
Rather  than  choosing  either  AE  or  <f>  to  be  the  inde¬ 
pendent  variable,  it  is  better  to  find  the  curve  which 
minimizes  the  total  square-error  in  both  the  vertical 
and  horizontal  directions  simultaneously.  We  have 
used  the  “simplex”  method  for  this  purpose  (specifi¬ 
cally.  a  variation  of  the  Amoeba  subroutine  by  Press 
rt  al.  (1986)),  since  it  can  minimize  by  trial-and- 
error  any  arbitrary  function. 


The  results  of  fitting  three  different  equations  to 
the  data  are  shown  in  Figure  11.  I'he  long-dashed 
lines  on  the  graphs  show  the  result  of  fitting  the  data 
to  a  straight  line: 

A  E  a  -t  f>4>  .  (16) 

But  the  data  appear  to  have  a  non-linear  trend.  T  his 
can  be  expected,  since  the  electrojet  current  depend:, 
on  both  the  potential  drop  and  the  Hall  conductiv 
ity,  and  precipitation-induced  conductivity  enhance¬ 
ments  are  related  to  the  potential  drop.  Therefore, 
it  is  not  unreasonable  to  test  relationships  with  an 
exponent  other  than  one.  The  other  equations  which 
have  been  fit  to  the  data  are:  a  power  law. 

AE  a  ■+  f>< t>r  .  (17) 

which  are  shown  with  the  solid  lines;  and  a  second- 
order  polynomial, 

AE  a  +  b<t>  4  c<J>‘  .  (IS) 

which  are  shown  with  the  short -dashed  lines.  I  he 
values  of  the  constants  which  were  determined  bv 
these  fits  are  given  in  Table  2.  The  total  square- 
error,  \2.  and  the  standard  deviation,  <?,  are  also 
included  in  Table  2.  The  quantity  \2  is  minimized 
by  the  curve-fitting  routine.  It  is  defined  for  our 
purposes  as: 
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FIG  11  SCATTER  PLOTS  OF  AE  INDEX  VS.  POLAR  CAP  POTENTIAL  AND  THE  RESULTS  OF 
FI  I  I  INC  CU  RV  ES  TO  THE  DATA.  These  data  include  only  the  satellite  passes  at  MLT  less  than  8  hours 
in  the  dawn  auroral  zone  and  less  than  20  hours  in  the  dusk  aurora)  zone,  (a)  Winter  period  in  the  northern 
hemisphere  (b)  Summer  period.  The  long-dashed  lines  show  the  fits  of  Equation  16,  a  straight  line.  The 
solid  lines  show  the  fits  of  Equation  17,  a  power  law.  The  short-dashed  lines  show  the  fits  of  Equation  IS.  a 
second-order  polynomial.  1  he  parameters  for  these  lines  are  given  in  Table  2. 
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(19) 


Note  that  the  AK  and  4*  errors  have  been  normal¬ 
ized  to  dimensionless  quantities  by  dividing  by  their 
maximum  values.  We  have  used  1400  nT  and  180 
k\\  the  maximum  values  on  the  plot  scales.  The 
standard  deviation  is  calculated  as: 


<T 


(20) 


These  quantities  have  some  similarities  to  the  con¬ 
ventional  definitions,  but  there  are  differences. 

The  power-law  formula  has  the  best  fit  to  the  data 
in  Figure  I  la,  the  winter  data,  and  the  second-order 
polynomial  has  the  lowest  error  in  Figure  11b,  the 
summer  data.  However,  there  is  not  a  very  large 
difference  in  the  errors.  The  standard  deviations  for 
the  non-linear  fits  are  not  significantly  different  from 
those  for  the  straight  line,  so  the  straight  lines  seem 
to  be  sufficient  for  modeling  the  data. 

By  inverting  Equation  16  and  using  the  values 
given  in  Table  2  for  a  and  b,  it  is  found  that 
♦  (kV)  26.8  4  0.152AE(nT)  in  the  winter  and 
<fr(kV)  -  19.2  4  0.116AE(nT)  in  the  summer.  These 
results  are  similar  to  those  found  by  Reiff  el  al. 
(1981)  (4>(kV)  41  +  0.1 1  AE(nT)),  Ahn  rt  al. 

(1984)  (<J>(kV)  =  36  4-  0.089AE(nT)),  and  by  Rich¬ 
mond  el  al.  (1990)  (<t>(kV)  22  +  0.1 19A E(nT)). 

One  problem  with  these  equations  is  that  they  give 
nonsensical  results  for  AE  at  very  low  potentials. 
The  earlier  studies  had  very  little  data  with  polar 
cap  potentials  below  30  kV,  and  our  own  data  show 
that  AE  is  very  low  when  the  potential  is  below  this 
value.  This  is  not  too  surprising,  since  the  auro¬ 
ral  zones  may  be  contracted  to  a  position  northward 
of  the  magnetometer  stations  when  the  convection  is 
weak.  Additionally,  the  low  potentials  will  usually  be 


associated  with  northward  IMF  11 7  conditions,  with 
erratic  convection  patterns  (Heppner  and  Maynard, 
1987)  and  disorganized  elertrojet  currents. 

Considering  that  Ahn  el  al.  (  1984)  and  Richmond 
el  al.  (1990)  ..ad  used  ground-based  measurements 
and  computer  algorithms  to  calculate  the  polar  cap 
potential  drops,  it,  is  remarkable  that  their  results  are 
in  such  close  agreement  with  these  obtained  bv  in 
situ  satellite  measurements.  This  indicates  that  the 
Kamide  method  is  fairly  accurate.  In  another  paper 
Ahn  el  al.  (1989)  had  also  used  the  magnetogram 
inversion  technique,  and  where  both  AE  indices  and 
calculated  polar  rap  potentials  were  given  their  data 
points  from  the  month  of  July  fit  within  the  spread  of 
our  points  in  Figure  7b.  They  had  points  at  (81  kV. 
562  nT)  and  (66  kV.  142  nT)  before  substorms  and 
at  (81  kV.  1062  nT)  and  (109  kV.  800  nT)  just  after 
substorm  maximum  phase  (the  last  AE  value  was  es¬ 
timated  from  a  graph  rather  than  given  in  the  text  ). 
Since  these  points  have  a  spread  similar  to  ours,  it 
appears  that  the  scatter  in  the  plots  is  more  an  in¬ 
herent  property  of  the  potential- AE  relation  than  a 
problem  with  our  satellite  sampling. 

Figures  5,  6,  and  7  have  trends  which  indicate 
that  the  ionospheric  conductivity  increases  as  the 
polar  cap  potential  increases,  since  the  AE  indices 
increase  more  rapidly  (i.e.,  the  slope  is  greater)  as 
the  potential  increases.  This  enhanced  conductivity 
is  due  to  electron  precipitation  associated  with  field- 
aligned  currents.  The  change  in  the  slope  is  much 
more  noticeable  in  AL  than  it  is  in  Al!,  which  is 
consistent  with  the  conductivity  in  the  dawn  auroral 
zone  increasing  much  more  than  in  the  dusk  auro¬ 
ral  zone  as  magnetic  activity  increases  (Hardy  el  al.. 
1987).  The  conductivity  enhancements  cause  an  ap¬ 
parent  discontinuity  in  the  A  E-potential  relations  in 
Figure  7,  occurring  at  about  69  kV  in  the  winter 
and  40  kV  in  the  summer.  This  difference  between 
the  seasons  is  curious.  One  possible  explanation  is 


TABLE  2.  RESULTS  OF  THE  CURVE  FITS  SHOWN  ON  FIGURES  1 1 A  (winter)  AND  IIP  (summer) 
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Equation 
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16 

176 

6.56 

3.942 

0.0967 

winter 

17 

65. 

0.79 

1.40 

3.255 

0.0878 

18 

8. 

0.45 

0.044 

3.351 

0.0891 

16 

165. 

8.60 

— 

3.536 

0.1  171 

17 

69. 

0.81 

1 .48 

3.550 

0.1173 

18 

67. 

5.13 

0.024 

3.447 

0.1 156 

summer 
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that,  due  to  tlie  greater  solar-induced  conductivity 
in  the  summer,  for  a  given  potential  t he  Pedersen 
and  field-aligned  currents  in  the  summer  hemisphere 
are  greater  than  in  the  winter  hemisphere.  There¬ 
fore,  at  a  given  potential  the  particle  precipitation 
and  the  precipitation-induced  conductivity  will  also 
be  greater  in  the  summer  hemisphere. 

Much  of  the  scatter  in  these  data  is  due  to  the  dy¬ 
namics  of  magnetospheric  substorms.  In  fact,  since 
substorms  may  be  due  to  an  “unloading'  of  stored 
energy  from  the  magnetotail  (Akasofu,  1989)  there 
may  be  times  when  substorms  occur  while  the  po¬ 
lar  cap  potential  (which  is  controlled  by  the  IMF) 
is  nearly  constant.  Therefore  we  should  expect  that 
at  some  potentials  the  AE  indices  will  have  a  wide 
range,  due  to  enhanced  electrojet  currents  during 
substorms.  The  role  of  substorms  is  evident  in  Fig¬ 
ures  9  and  10.  Wherever  a  substorm  growth  phase 
has  been  identified  the  AE  indices  are  down  near  the 
lower  10%  curve:  during  substorm  expansion  phase 
the  AE  indices  move  upward  toward  the  90%  curve; 
and  during  the  recovery  phase  the  AE  indices,  and 
sometimes  the  polar  cap  potential,  start  to  move 
downward. 

Kan  ft  al.  (1988)  have  proposed  that  the  polar  cap 
potential  must  reach  a  certain  critical  level  (70  kV 
in  their  model)  in  order  for  substorm  onset  to  occur. 
Our  results  tend  to  support  that  aspect  of  their  the¬ 
ory.  Figure  7  shows  that  below  approximately  00  kV 
the  AE  indices  tend  to  remain  low.  with  little  verti¬ 
cal  scatter,  but  above  60  kV  the  AE  indices  become 
large  and  widely  scattered,  indicative  of  substorm 
activity.  Figures  9  and  10  are  even  more  convincing: 
the  potentials  during  the  growth  phase  are  between 
10  and  80  kV  and  the  potentials  during  substorm  ex¬ 
pansion  are  above  60  kV.  During  substorm  recovery 
several  of  the  measured  potentials  are  below  80  kV, 
but  there  are  also  instances  where  the  potentials  are 
still  high.  According  to  the  Kan  cl  al.  (1988)  model, 
substorm  recovery  occurs  when  either  the  polar  cap 
potential  drops  below  the  critical  level  or  the  diffuse 
auroral  conductance  belt  moves  equatorward  from 
the  convection  reversal. 


4.  CONCLUSION 

We  have  shown  that  the  AU,  AL,  and  AF,  indices 
can  be  related  to  the  respective  dusk  auroral,  dawn 
auroral,  and  polar  cap  potentials.  Due  to  the  large 
number  of  variables  the  relationships  are  approxi¬ 
mate  rather  than  exact.  Linear  relationships  which 
indicate  the  bounds  of  80%  of  the  data  points  and 
median  slopes  are  specified.  Different  parameters 
are  required  in  order  to  bracket  the  data  at  different 


times  of  the  year.  Since  the  AE  indices  are  measured 
only  In  one  hemisphere,  these  indices  are  greater  dur¬ 
ing  the  local  summer,  when  the  ionospheric  conduc¬ 
tivity  is  greater  than  in  the  winter. 

At  very  low  polar  cap  potentials  the  API  measure 
merits  do  not  appear  to  be  very  accurate,  due  to  the 
movement  of  the  auroral  oval  northward,  away  from 
the  magnetometer  stations.  At  high  potentials  the 
electrojet  currents  are  influenced  bv  precipitation- 
induced  conductivity  enhancements.  The  enhanced 
conductivities  are  more  noticeable  in  the  AL  index 
than  in  the  AU  index,  since  the  westward  electrojet 
is  associated  with  upward  current  and  precipitating 
electrons  while  the  eastward  electrojet  is  associated 
with  downward  current. 

Our  data  are  consistent  with  the  prior  observa 
tions  of  Katnide  and  Baumjohann  (1985),  Ahn  rl  al. 
( 1 989),  and  many  others  that  t he  increase  in  t he  eler- 
trojet  current  intensity  prior  to  substorm  onsets  is 
due  to  an  increase  in  the  electric  field/potential.  and 
conductivity  enhancements  after  onsets  contribute  to 
the  substorm  cnTents.  We  also  find  that  substorms 
appear  to  occur  when  the  polar  cap  potent  ial  exceeds 
a  threshold  of  approximately  60  kY.  which  supports 
(but  does  not  necessarily  prove)  a  theory  by  Kan  il 
al  (1988)  for  the  cause  of  substorm  onsets. 

Although  the  AE-potential  graphs  do  have  a  lot 
of  “scatter",  there  is  sufficient  order  to  the  data 
to  derive  a  least  squares  fit  equation.  By  invert¬ 
ing  the  equation  it  is  possible  to  derive  reasonable 
estimates  for  the  potentials  from  the  AF.  indices. 
This  would  be  a  logical  option  whenever  an  esti¬ 
mate  of  the  polar  cap  potential  is  required  for  a 
specifir  time  period  for  which  AE  is  known.  Al¬ 
though  the  AE  indices  do  have  faults  (Akasofu,  1989) 
they  are  more  accurate  than  Kp  for  indicating  the 
state  of  the  magnetosphere-ionosphere  system  at  a 
given  moment.  Future  modelers  of  ionospheric  con 
vection,  conductivities,  currents,  precipitation,  and 
Joule  heat ing  should  index  t heir  models  to  A E  rather 
Ilian  Kp. 

Seasonal  effects  also  need  to  be  taken  into  consid¬ 
eration  more  frequently  than  they  have  in  the  past. 
In  other  words,  most  people  who  have  used  the  A  F. 
imiex  in  one  way  or  another  have  done  so  with  little 
regard  for  the  time  of  the  year,  it  could  be  useful 
to  have  vet  another  index,  similar  to  AE,  but  with 
adjustments  made  for  annual  variations. 
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